Incompressible liquid state of rapidly-rotating bosons at filling factor 3/2 



in 
o 
o 

(N 



E.H. Rezayi 1 , N. Read 2 and N.R. Cooper 3 
1 Department of Physics, California State University, Los Angeles, CA 90032 
2 Department of Physics, Yale University, P.O. Box 208120, New Haven, CT 06520-8120 and 
3 T.C.M. Group, Cavendish Laboratory, Madingley Road, Cambridge, CBS OHE, United Kingdom 

Bosons in the lowest Landau level, such as rapidly-rotating cold trapped atoms, are investigated 
numerically in the specially interesting case in which the filling factor (ratio of particle number to 
vortex number) is 3/2. When a moderate amount of a longer-range (e.g. dipolar) interaction is 
included, we find clear evidence that the ground state is in a phase constructed earlier by two of us, 
in which excitations possess non-Abelian statistics. 
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There is increasing interest in rapidly-rotating ultra- 
cold atoms in a trap. For bosons rotating at a moderate 
frequency, a vortex lattice is observed Q . At sufficiently 
high rotation frequency (close to the natural frequency of 
a harmonic trap), it is expected that the lattice melts and 
is replaced by a series of highly-correlated liquids that can 
be related to fractional quantum Hall (QH) states. The 
crucial parameter is the ratio v = N/Ny of the number 
of atoms, N, to the number Ny of quantized vortices 
that would pierce the cloud if it were a Bose condensate. 
This ratio corresponds to the Landau level filling factor v 
in the fractional quantum Hall effect. Previous theoreti- 
cal work on this regime has emphasized the importance 
of the lowest Landau level (LLL) when interactions are 
weak and the temperature is low, and pointed out that in 
this restricted space of states the Laughlin v = 1/2 state 
is the exact ground state for the standard "contact" 
form of interaction Q. Later, evidence of a sequence of 
correlated liquids was found at v = fc/2, k = 1, 2, 3, 
. . . for v < Vpjm, and a vortex lattice at v > v c , with 
v c ~ 6-10 0^ U The correlated liquids were found Q 
to have large overlaps with states proposed by two of us 
(RR) some time ago 0. The k = 1 liquid is the famil- 
iar Laughlin state, while k = 2 is the Moore-Read (MR) 
paired state Q. More recent work strengthens the case 
for the MR state at v = 1, and provides some evidence 
for liquid states at still other filling factors not in the 
sequence v = k/2 §|, however results at v = 3/2 were in- 
conclusive. Meanwhile, experiments are approaching the 
LLL regime 9], though the filling factors are still 3> v c . 
Very recently, Bose condensation has been achieved in 
atoms with a large magnetic dipole moment [13, and 
the effect of dipolar interactions on the vortex lattices 
and on the quantum fluids at low filling factors has been 
investigated [ill ]. 

In this paper, we study numerically the next mem- 
ber, v = 3/2, of the RR sequence for system sizes larger 
than in Refs. jj, la We consider the s-wave (con- 
tact) interaction in the LLL, and the effect of adding 
a longer-range component such as the dipolar interac- 
tion (for dipole moments oriented parallel to the rota- 
tion axis). This case, k — 3, is of interest for sev- 
eral reasons: For k > 1, the quasiparticle excitations 



of each of the sequence of states introduced in Ref. 
have the fascinating property of non-Abelian statistics 
7]. This makes these states even more exotic than the 
Laughlin and hierarchy/composite-fermion ^1 states, 
in all of which the quasiparticles have fractional, but 
Abelian, statistics, k — 3 is the smallest k value in the 
RR sequence for which the non-Abelian statistics sup- 
port universal quantum computation |l3j |. For this fill- 
ing factor, v = 3/2, there is also an alternative candi- 
date, which is a hierarchy/composite-fermion phase. Us- 
ing both the sphere and the torus geometries (to avoid 
edge effects), we find clear evidence that the boson sys- 
tem at v = 3/2 with a moderate amount of longer-range 
interaction added is in the phase described by our trial 
states, and hence is non-Abelian; this may also be the 
case for the pure s-wave interaction, as suggested by Ref. 
. The evidence comes from the energy spectrum on the 
torus, which shows a ground-state doublet with very large 
overlaps with the RR trial states, and a relatively large 
gap to higher excited states, and from the two-particle 
correlation function. 

The conventional effective interaction Hamiltonian for 
the atoms, representing s-wave scattering at low momen- 
tum, is 
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where g = AnT^a/M (a is the s-wave scattering length, 
and M is the mass of an atom). It is of interest to con- 
sider also electric or magnetic dipole interactions between 
the atoms, of the form 
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where = (r, — i*j)/|rj — r 3 -|, and the p^s are unit vec- 
tors representing the (fixed) dipole moments. We will 
assume that the dipole moments are parallel to the 3 
axis; then any (5-function term that may accompany the 
dipolar interaction can be absorbed into the s-wave in- 
teraction term in H m t = H s + Hd\ P ■ 

We work in an axially symmetric harmonic trap, with 
frequencies U3, uj±_ for motion respectively along, and 
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perpendicular to, the symmetry (3-) axis. For a system 
of bosons rotating rapidly about the 3-axis, we will re- 
strict the atoms to the lowest Landau level (LLL) states, 
which are the single-particle states of lowest energy for 
each value to > of angular momentum L^/h (and thus 
have no excitation along the 3-axis). This is based on 
the assumption that interactions are weak, that is gn and 
Can (n is the typical density of particles in the drop of 
atoms) are small compared with the energy for excitation 
to higher states, that is with 2fiwj_ and Huj^ HH. The 
LLL states can be represented by functions in two dimen- 
sions, u m (z) — z m e~\ z \ / 4 /\/27r2 m TO!, where z — x + iy 
(henceforth we use units in which the magnetic length £g 
and h are 1 0)- 

When working on the sphere 0], there is a finite 
number Ny of flux quanta penetrating the surface of ra- 
dius R = and there are Ny + 1 single-particle 
states, which have wavefunctions u m (z) cx z m /(l + 
\z\ 2 /AR 2 ) 1+Nv / 2 , m = 0, N v in terms of the co- 
ordinate z in the plane (from stereographic projection of 
the sphere). In the torus geometry (i.e. periodic bound- 
ary conditions on a parallelogram), there are Ny LLL 
basis states when Ny flux pierce the system. 

In the subspace of LLL states in the infinite plane, the 
interaction Hamiltonian H{ nt can be represented by the 
pseudopotentials V m , to = 0, 2, ... . V m is the interac- 
tion energy for a single pair of bosons of relative angular 
momentum to (only even to are relevant for bosons) |l5| . 
The pseudopotentials can be obtained 0] from the ana- 
lytic expressions in the limit (for simplicity) of vanishing 
thickness of the two-dimensional fluid, l%/l±_ — ► In 
this limit an infinite constant has to be absorbed into g. 
The to 7^ pseudopotentials are determined entirely by 
the dipolar interaction, and hence their ratios are fixed; 
they are plotted as an inset in Fig. [3 below. Vq can be 
treated as an independent parameter, so that the ratio 
of V2/V0 is a dimensionless parameter characterizing the 
interaction apart from one overall energy scale. We also 
studied an interaction in which V m for to > 2 is set to 
zero. This description of the interaction can be extended 
to the sphere (also using rotation symmetry), and to the 
torus. 

A p-body (5-function interaction, 
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projected to the LLL is also of interest (this form is cor- 
rect for the plane and torus geometries, and for the sphere 
there is a corresponding rotationally-invariant form). 
The parafermion states found in ref. are unique, exact 
zero-energy eigenstates of such interactions for p = k + 1 
when N is divisible by k and Ny = 2N/k — 2 (on the 
sphere), so that v = lhxiTv^oo N /Ny = k/2. These states 
serve as trial wavefunctions with which the exact ground 
states for H nrt can be compared. They represent incom- 
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FIG. 1: Low-lying spectrum for 18 bosons for dipolar model 
Hi n t with V2/V0 — 0.380 on the torus vs. the pseudomomen- 
tum K — IKI. for two different unit cells. 



pressible liquid phases, in which the excitations enjoy 
non-Abelian statistics for k > 1 (however, we will later 
suggest that a caveat to this statement is required). 

Before turning to our results, we review some as- 
pects of the LLL on a torus |l6(. On the torus at 
v = N/Ny = 3/2, translational symmetry implies that 
all energy eigenstates possess a trivial center-of-mass de- 
generacy of 2, which is exact for any size system, and 
also that there is a conserved pseudomomcntum K 
which is a vector lying in a certain Brillouin zone. In the 
RR phases, the ground states have a net degeneracy k + 1 
in the thermodynamic limit, which is connected with the 
non-Abelian statistics 0,0 ■ For k = 3, this 4-fold degen- 
eracy is made up of the trivial factor 2 (which is always 
discarded in numerical studies), together with a further 
2-fold degeneracy, which in general becomes exact only 
in the thermodynamic limit; all these ground states have 
K = 0. (For the 4-body interaction, the fourfold degen- 
eracy is exact for any size, as the ground states have ex- 
actly zero energy.) By contrast, the hierarchy/composite- 
fermion ground state for bosons at v = 3/2 possesses only 
the trivial 2-fold degeneracy. Then for an incompressible 
fluid on the torus, the spectrum of a sufficiently large 
system in one of these two phases should exhibit a nearly 
degenerate pair of ground-states (resp., a single ground 
state) at K = 0, separated by a clear gap from a region 
of many states at higher energy eigenvalues. 

In Fig. ^ we show the spectrum for 18 particles on 
the torus for H int with V2/V0 — 0.380 (energies are in 
units of Cd for the dipolar H). A clear doublet can be 
seen at K = 0, and the rest of the spectrum is well- 
separated compared with either the spacing of the dou- 
blet, or the spacing among the levels above the gap. This 
spectrum leaves little doubt that the fluid is incompress- 
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FIG. 2: The total overlap- squared as a function of V2/V0 
on the torus (square unit cell), for 18 particles, for both the 
dipolar and the Vq, V-j-only model interactions. The inset 
shows the m/0 pseudopotentials for the dipolar interaction. 



FIG. 3: Low- lying spectrum for 21 particles on the sphere 
for dipolar interaction with V2/V0 ~ 0.350 vs. equivalent 
wavevector K — L/R (L is the angular momentum). Inset 
shows the overlap- squared, with symbols as in Fig. |5] 



ible. The results are similar for the two geometries of the 
torus (hexagonal and square unit cells) shown. Spectra 
at nearby V2/V0, and at other aspect ratios, are similar. 

The overlap of each of the two low-lying K = states 
with the trial ground state doublet of the 4-body inter- 
action can be calculated. We may add the squares of 
these overlaps to obtain the total overlap-squared (which 
is at most 2) of the two-dimensional subspace of i?i n t 
with that of H4. This total overlap-squared is plotted 
versus V2/V0 for both the Vb plus dipole interaction and 
the Vq, V2-only interaction in Fig. (the total is made 
up of roughly equal contributions from each of the two 
low- lying states, across the whole range of V2/V0 val- 
ues). For sufficiently large V2/V0, both cases show a very 
abrupt drop in overlap which is due to a phase transi- 
tion. Beyond this point, the low- lying spectrum shows 
indications that the ground state in the thermodynamic 
limit breaks translational symmetry. The transition is 
probably first-order, but does not occur here by crossing 
of energy levels of different symmetry, because a crystal 
and the RR fluid both have ground states at K = 0. This 
transition is similar to that found in other recent work 
at v = 1/2 [13 and will not be pursued here. We note, 
however, the similarity of these overlap curves with those 
for v = 1/2 in the case of fermions as the interaction is 
changed, which show similar trends and similar curves 

The N — 18 overlaps should be viewed as significantly 
large. Note that for a random vector in D dimensions, the 
probability of obtaining an overlap with a given vector 
(or two-dimensional subspace) greater than some given 
value is of order e~ D as D — > 00. The probability of 
obtaining similar overlaps for two random vectors is of 
order the square of this. In our case, the K = block of 



the matrix contains (with reflection symmetry included) 
about 242,000 states. However, there are other symme- 
tries that we did not use, such as point operations, that 
commute with translations at the zone center and would 
further reduce somewhat the dimension of the relevant 
space. The large value of not only the ground state, but 
also the lowest excited state, overlap with the trial sub- 
space is very strong evidence that these systems are in 
the RR phase, and would not support an interpretation 
as hierarchy/composite-fermion states. 

Fig. |21 shows the spectrum in the spherical geometry, 
with V2/V0 = 0.350. Again, a clear gap separates the 
ground state from the rest of the spectrum. The results 
are consistent with those on the torus. 

The two-particle correlation function g(r) for the 
ground state on the sphere is shown in Fig. for the Vq , 
V2-only interaction for 21 particles, together with that 
of the RR trial state (the ground state of H4) for 18 
particles, for comparison. In general, g(r) in a strongly- 
correlated state should show a correlation hole at short 
distances (i.e. g(0) < 1), and in an incompressible fluid 
state it should tend to 1 exponentially in r at long dis- 
tances. In the present case, while g(r) is less than 1 at 
r = 0, it has a local maximum there (mentioned in Ref. 
@ , and this seems to be a general feature of the ground 
states for v > 1/2). This surprising result may indi- 
cate a general tendency for bosons to cluster, at least in 
these rotationally-invariant ground states. The behavior 
at large r suggests that these systems are on the verge 
of being larger than the relevant correlation length. For 
the RR trial state, g(r) for 18 particles has essentially 
converged to its thermodynamic limit. 

We can gain insight into the structure of g{r) as v 
increases by considering the vortex lattice. In this state 
there are holes in the density at some fixed positions. Ac- 
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FIG. 4: Two-particle correlation function g(r) of various 
ground states on the sphere, plotted against great-circle dis- 
tance r. N = 21 g(r) is for a V0-V2 model. 



FIG. 5: Low-lying spectrum (in units of g) for 18 particles 
on the torus vs. pseudomomentum K, for H a for two different 
unit cells. 



cordingly, g(r) in the vortex lattice, averaged over shifts 
and rotations of the lattice to represent the finite-size sys- 
tem, will have maxima and minima at arbitrarily large 
r; in particular, there will be a maximum at r = 0. It 
seems that in the RR series, the oscillations in g{r) die 
out with increasing r at least for small k, but that as 
k increases the oscillations extend out to larger r, and 
the maximum at r = becomes larger. We suspect that 
for sufficiently large k, the RR trial states actually ex- 
hibit a transition to vortex-lattice long-range order (or 
possibly to one of the other ordered states found in Ref. 

This would be analogous to the Laughlin states 
at v = 1/m, which exhibit increasing oscillations with 
increasing to, and eventually develop crystalline order of 
the particles (g(0) = in these cases) |2(. Note that in 
both cases, the long-range order in the trial states is (or 
would be) occurring for the exact ground states of some 
family of special Hamiltonians, that are not the ones of 
most physical interest, yet the transition is of the same 
type as that for the latter (though the critical v may be 
much different, as for the Laughlin states). 

Fig. shows the spectrum of the pure contact interac- 
tion H s (i.e. V2/V0 = 0) on the torus. This, in conjunc- 
tion with the corresponding g(r) on the sphere for N = 24 
particles shown in Fig. shows no clear signal that the 
system is incompressible. While it is possible that the 
correlation length is simply larger, and that larger sys- 
tems would exhibit incompressibility, we cannot rule out 
the possibility that the fluid is this region is compressible, 
or that there is a breaking of translational symmetry in 
the thermodynamic limit. Indeed, the dependence of the 
spectrum on the aspect ratio in the torus geometry sug- 
gests the presence of a broken symmetry "stripe" phase 
competing with the RR phase . 

In conclusion, we have exhibited clear evidence that 
for bosons in the lowest Landau level at filling factor 



3/2, when a moderate amount of longer-range interaction 
is included, the ground state is an incompressible fluid 
of a type that possesses non-Abelian statistics for the 
quasiparticle excitations. 
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